Deconfined criticality for the S = 1 spin model on the spatially 

anisotropic triangular lattice 



Yoshihiro Nishiyama 

Department of Physics, Faculty of Science, 
Okayama University, Okayama 700-8530, Japan 
(Dated: January 12, 2011) 

Abstract 

The quantum S = 1 spin model on the spatially anisotropic triangular lattice is investigated 
numerically. The nematic and valence-bond-solid (VBS) phases are realized by adjusting the 
spatial anisotropy and the biquadratic interaction. The phase transition between the nematic and 
VBS phases is expected to be a continuous one with unconventional critical indices (deconfined 
criticality). The geometrical character (spatial anisotropy) is taken into account by imposing the 
screw-boundary condition (Novotny's method). Diagonalizing the finite-size cluster with < 20 
spins, we observe a clear indication of continuous phase transition. The correlation-length critical 
exponent is estimated as = 0.92(10). 
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I. INTRODUCTION 



According to the deconfined-criticality scenario,-"- in two dimensions, the phase transi- 
tion separating the valence-bond-sohd (VBS) and antiferromagnetic phases is continuous; 
naively,^ such a transition should be discontinuous, because the adjacent phases possess 
distinctive order parameters such as the VBS coverage pattern and the sublattice magneti- 



zation, respectively. A good deal of theoretical 
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investigations has been made to support 
it was claimed that the transition would be 



this scenario. (On the contrary, in Refs. 
a weak first-order one.) 

The magnetic frustration is a clue to the realization of the VBS phase. Actually, the 
square-lattice antiferromagnet with the next-nearest-neighbor interaction ( J1-J2 model) ex- 
hibits the VBS phase around the fully frustrated {J2/ Ji ~ 0.5) regime.— The quantum 
Monte Carlo method is not applicable to this problem because of the negative-sign prob- 
lem. So far, the J1-J2 model has been studied with the series-expansion^"^ and numerical- 
diagonalization^ methods. 

Alternatively, one is able to realize the VBS phase through incorporating the biquadratic 
interaction^^i^; correspondingly, one has to enlarge the magnitude of spin to S > 1/2. 
The biquadratic interaction (unlike the magnetic frustration) is tractable with the quantum 
Monte Carlo method. The deconfined criticality is realized by tuning the spatial anisotropy.— 
That is, as the interchain interaction increases, a transition from the VBS phase to either ne- 
matic or antiferromagnetic phase occurs.— Meanwhile, it turned out that the ring-exchange 
(plaquette-four-spin) interaction also induces the VBS phase even for S = 1/2. Extensive 
Monte Carlo simulations^!^ support the deconfined-criticality scenario; the results are com- 
pared with ours in Sec. IIVI The transition occurs for a considerably large ring exchange; the 
antiferromagnetic order would be far more robust than that of VBS. As mentioned above, 
the character of the singularity is controversial^^—; possibly, the log corrections^^ affect the 
scaling analysis. Even for the S = 1/2 spatially anisotropic triangular antiferromagnet,— 
The ring exchange gives rise to the VBS phase.— This model is intractable with the quan- 
tum Monte Carlo method, and analytical considerations provide valuable information as to 
the deconfined criticality.—"— 

In this paper, we investigate the S = 1 spatially-anisotropic-triangular-lattice model 
with the biquadratic interaction by means of the numerical diagonalization method. To be 



specific, the Hamiltonian is given by 

H = -JEb's. ■ s, + (s. ■ s,f] -j'Y. (s. ■ s,)l (1) 

fo-> m) 

Here, the quantum S' = 1 spins {Sj} are placed at each triangular-lattice point i] see Fig. [T] 
(a). The summation runs over all possible nearest-neighbor (skew-diagonal) 

pairs. The parameter J {J') denotes the corresponding coupling constant. Hereafter, we 
consider J' as the unit of energy (J' = 1). Along the J bond, both quadratic and biquadratic 
interactions exist, and the parameter j controls a strength of the former component. The J'- 
bond interaction is purely biquadratic. The interaction J interpolates the one-dimensional 
(J = 0) and square-lattice {J' — t- oo) structures. Correspondingly, the VBS and spin-nematic 
phases appear, as the interaction J varies; see a schematic phase diagram. Fig. [2l In order to 
take into account such a geometrical character, we implement the screw-boundary condition 
[Fig. [11(b)] through resorting to Novotny's method (Sec. HT]) . 

The rest of this paper is organized as follows. In Sec. [TTl we explain the simulation 
scheme. We also make an overview on the biquadratic-interaction spin models relevant to 
ours. In Sec. IIIIl we demonstrate that the present model exhibits a clear indication of 
deconfined criticality at a moderate value of J. We also analyze the criticality with the 
finite-size-scaling theory. In Sec. llVt we present the summary and discussions. 

II. SCREW-BOUNDARY CONDITION: NOVOTNY'S METHOD 

In this section, we present the simulation scheme (Novotny's method^i^). A brief 
overview on the biquadratic-interaction spin models follows. 

Before commencing a explanation of technical details, we present a basic idea of Novotny's 
method. We implement the screw-boundary condition for a finite cluster with spins; see 
Fig. [11(b). Basically, the spins, {Sj} {i < N), constitute a one-dimensional {d = 1) structure, 
and the dimensionality is lifted to = 2 by the bridges over the long-range pairs. The present 
system ([T|) has a spatial anisotropy governed by J. We take into account such a geometrical 
character through imposing the screw-boundary condition. According to Novotny, the long- 
range interactions are introduced systematically by the use of the translation operator P; 
see Eq. ([3]), for instance. The operator P satisfies the formula 

-Pi'S'i, 5*2, . . . , Sn) = \Sn, 5*1, . . . , ^jv-i). (2) 
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Here, the base |{5'i}) diagonalizes each of {S^}; namely, the relation SI.\{Si}) = Sk\{Si}) 
holds. 

Novotny's method was adapted to the quantum 5 = 1 XY model in c? = 2 dimensions.— 
Our simulation scheme is based on this formalism. In the following, we present the modifi- 
cations explicitly for the sake of selfconsistency. The XY interaction Hxy, Eq- (4) of Ref. 
, has to be replaced with the Heisenberg interaction 
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N 



Hxxx{v) = ^(P^Sfp-'^Sf + P^'Sfp-^'Sf + P^'SlP-^'St). (3) 



i=l 



Additionally, we introduce the biquadratic interaction 

1 I N 5 

H,{v) = --Hxxxiv) + - E E P''Q?P-''Qt- (4) 

^ ^ i=l a=l 

The definition of {Q"} and an algebra are presented in the Appendix. Based on these 
expressions, we replace Eq. (3) of Ref. 
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with 



n = -JljHxxxiVN) + jHxxxi^ - 1) + H,{^m) + H,i^/N - 1)] - J'H.il). (5) 

We diagonalize this matrix for < 20 spins in Sec. Illli The above formulae complete the 
formal basis of our simulation scheme. However, in order to evaluate the above Hamiltonian- 
matrix elements efficiently, one may refer to a number of techniques addressed in Refs. 
and 36. 
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Last, we overview the biquadratic-interaction spin models. As mentioned in the Intro- 
duction, the present model reduces to the one- dimensional and square-lattice models in the 
limiting cases J = and J — t- oo, respectively. Each of these limiting cases has been 
studied extensively. Here, we devote ourselves to the nearest-neighbor interaction of the 
generic form, cos6'Si ■ S-,- + sin6'(Sj ■ S-,)^, parameterized by 6. The regime tt < 6 < 27r 
is relevant to the present research. As for d = 1 (Ref. 1371), the dimer phase appears in 
57r/4 < 6 < Ttt/A. Namely, around 6 ^ 37r/2, the stability of the dimer (VBS) phase would 



be maximal. For d = 2 (Ref. l38l ). the ferromagnetic, nematic, and antiferromagnetic phases 
appear in 6 < 5tt/A, dn/A < 9 < 37r/2, and 37r/2 < 9, respectively. The phase diagram in 
Fig. [2] is based on these preceding studies. 



III. NUMERICAL RESULTS 



In this section, we present the simulation results. We employed the simulation scheme 
developed in Sec. [TTl We devote ourselves to the analysis of the transition between the 
nematic and VBS phases (deconfined criticality) ; see Fig. |21 We treat a variety of system 
sizes = 10, 12, ... , 20 (A^ is the number of spins within a cluster). The linear dimension 
L of the cluster is given by 

L = Vn. (6) 

A. Critical point 

In this section, we investigate a location of the phase boundary separating the VBS and 
nematic phases. 

In Fig. 131 we plot the scaled energy gap LAE for various J and = 10, 12, . . . , 20. 
The quadratic-interaction strength j is fixed to j = 0.5. The symbol AE denotes the first 
excitation gap. According to the finite-size scaling, the scaled energy gap LAE should be 
scale-invariant at the critical point. In fact, we observe an intersection point at Jc ~ 0.28, 
which indicates an onset of the J-driven phase transition. As mentioned in the Introduction, 
the nature of this singularity is under current interest. The present result indicates that the 
singularity is a continuous one; the critical phenomena (below the upper critical dimension) 
should be described by the finite-size scaling. 

In Fig. m we plot the approximate transition point Jc{Li,L2) for [2/{Li + -^2)]^ with 
10 < A^^i < A'2 < 20 {Li^2 = \J Ni^2)- The parameters are the same as those of Fig. [3l Here, 
the approximate transition point denotes a scale-invariant point with respect to a pair of 
system sizes (Li,L2). Namely, the following relation holds; 

LrAE{L,)\j=j^(L,M) = L2AE{L2)\j=ML,M)- (7) 

The least-squares fit to the data of Fig. H] yields an estimate Jc = 0.285(5) in the thermo- 
dynamic limit L — )■ 00. In order to appreciate possible extrapolation errors, we made an 
alternative extrapolation with the 1/L^-abscissa scale. Thereby, we obtain Jc = 0.285(3). 
The extrapolation errors appear to be negligible. (As a matter of fact, we surveyed a wide 
range of j, and found that the parameter j = 0.5 yields an optimal finite-size behavior.) 
Hence, we obtain an estimate Jc = 0.285(5). 
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Making simular analyses for various values of j, we arrived at a phase diagram, Fig. |2j 
The phase boundary around j ~ and 1 is ambiguous because of finite-size errors. Possibly, 
around j ~ 1, the magnetic structure (antiferromagnetic order) conflicts with the screw- 
boundary condition, resulting in an enhancement of finite-size errors. On the one hand, in 
j < 0, the character of the transition changes to a discontinuous one, and the finite-size- 
scaling method becomes invalid. Hence, the finite-size behavior improves around the midst 
(j ^ 0.5) of < j < 1. 

Last, we address a number of remarks. First, in the scaling analysis. Fig. |3l we as- 
sumed the dynamical critical exponent z = 1, following the conclusion of the Monte Carlo 
simulations.— i2I Second, we argue a possible systematic error for the data Jc{Li, L2) in Fig. 
m As a matter of fact, for large system sizes, the data Jc{Li,L2) exhibit an enhancement, 
suggesting that the extrapolated value of Jc should be larger than 0.285(5). However, in the 
subsequent analyses, the extrapolated value of Jc is no longer used, and systematic devia- 
tions are less influential. Nevertheless, as to the singularity of Jc, it has to be mentioned 
that the present data cannot exclude a possibility of weak-first-order transition, for which 
the scaling approach becomes invalidated. 



B. Correlation-length critical exponent u 

In this section, we analyze the criticality found in Sec. IIII A[ 
In Fig. m we plot the approximate critical exponent 

J=Jc(Li,L2) 

for 2/(Li + L2) with 10 < Ni < N2 < 20. The parameters are the same as those of 
Fig. [31 The least-squares fit to these data yields u = 0.97(2) in the thermodynamic limit. 
This estimate may be affected by systematic errors. (The statistical error would be an 
underestimate.) In order to appreciate an error margin, we made an alternative extrapolation 
with the 1/L^-abscissa scale. Thereby, we obtain u = 0.86(1). The discrepancy indicates an 
amount of systematic errors. As a result, we arrive at an estimate 

i/ = 0.92(10), (9) 

which covers the above results obtained via independent extrapolations. 
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This is a good position to address a remark on the scahng behaviors of Figs. [3ll5l As 
mentioned in the Introduction, notorious log corrections were observed for the S = 1/2 
square-lattice antiferromagnet with the ring exchange. Our data, on the contrary, appear 
to exhibit moderate corrections to scaling, particularly, in Figs. E] and HJ For the former 
antiferromagnet, the VBS phase emerges for a considerably large ring exchange, indicating 
that the antiferro magnetic phase is robust. On the one hand, our model exhibits a stable 
VBS state owing to the spatial anisotropy (one-dimensionality in the J = limit), and 
the nematic phase turns into the VBS phase at a moderate coupling strength, ~ 0.3. 
We suspect that these peculiarities of the present (rather artificially designed) model bring 
about improved scaling behaviors. 

Last, we make a consideration of the abscissa scale of Fig. |5] (jlj). The critical exponent 
(point) has a leading correction of 0{L^^) [0(L^'^^^/'^)]; here, the symbol u denotes the 
index for corrections to scaling. At present, the index uj for the deconfined criticality is 
unclear. As a reference, one may refer to that of the = 3 Heisenberg universality class, 
oj = 0.773 (Ref. [sol. Making use of this value, we set the abscissa scale to that depicted in 
Fig. E]©. 



IV. SUMMARY AND DISCUSSIONS 

The S = 1 spin model on the spatially anisotropic triangular lattice, Eq. ([1]), was in- 
vestigated numerically. As the spatial anisotropy J changes, the VBS and nematic phases 
appear (Fig. [2]). Hence, this rather artificial model provides a candidate for the analysis of 
the deconfined criticality, which is arousing much attention recently. We employed Novotny's 
method (screw boundary condition) in order to take into account such a geometrical char- 
acter (spatial anisotropy). 

As a result, we observe a clear indication of the J-driven criticality through the finite-size- 
scaling analysis (Fig. [3]). This result supports the deconfined-criticality scenario. Thereby, 
we estimate the correlation-length critical exponent as u = 0.92(10). 

As a reference, we overview related studies. For the square-lattice antiferromagnet with 



the ring exchange, the estimates, u = 0.78(3) (Ref. l26|) and u = 0.68(4) (Ref. l27|), were 
reported. As for the spatially-anisotropic-triangular antiferromagnet with the ring exchange, 
the exponent z/ = 0.80(15) was obtained.— These results are to be compared with the index. 
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FIG. 1. (a) We consider a spatially anisotropic triangular lattice; the Hamiltonian is given by 
Eq. ([1]). The interaction J interpolate the one- and two-dimensional lattice structures in the 
limiting cases J = and J — )• oo, respectively, (b) In order to take into account such a geometrical 
character, we implement the screw-boundary condition. As shown in the drawing, a basic structure 
of the cluster is an alignment of spins {Sj} {i < N). Thereby, the dimensionality is lifted to d = 2 by 
the bridges over the (\/iV)-th neighbor pairs through the J bonds. Technical details are explicated 
in Sec. [Ill 



u = 0.7112(5) (Ref. |4y), for the d = 3 Heisenberg universality class. Our result indicates a 
tendency toward an enhancement for the correlation- length critical exponent, as compared 
to that for the d = 3 Heisenberg universality class. Taking the advantage of the numerical 
diagonalization method, we are able to extend the interactions so as to eliminate finite-size 
errors. A frustrated interaction along the J'-bond direction may stabilize (extend the regime 
of) the VBS phase substantially. This problem will be addressed in the future study. 



Appendix: A reduction formula for the biquadratic interaction 

The biquadratic interaction (Sj ■ Sj)^ reduces to a seemingly quadratic form 

(S, ■ S,)2 = -S, ■ S,/2 + ^ + 4/3. (A.l) 

a=l 
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FIG. 2. A schematic phase diagram for the 5 = 1 spatially-anisotropic-triangular-lattice model, 



37 and 



3^, 



Eq. (HI), is presented. The hmiting cases J — t- and J — t- oo were studied in Refs. 
respectively. The solid (dashed) lines stand for the phase boundaries of discontinuous (continuous) 
character. The dotted lines are ambiguous. We investigate the phase boundary separating the 
nematic and VBS phases. 




0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 



J 

FIG. 3. The scaled energy gap L/S.E is plotted for various J and N = 10, 12, ... , 20. The 
quadratic-interaction strength j is fixed to j = 0.5. (J' is the unit of energy.) We observe a clear 
indication of the deconfined criticality around J ~ 0.28. 
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FIG. 4. The approximate critical point Jc(Li, L2) ([ZD is plotted for [2/{Li + L2)f with 10 < iVi < 
A''2 < 20. The parameters are the same as those of Fig. [3l The least-squares fit to these data yields 
Jc = 0.285(5) in the thermodynamic limit. 
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0.35 



FIG. 5. The approximate critical exponent v{Li,L2) ([8]) is plotted for 2/(Li + L2) with 10 < 
A^i < N2 < 20. The parameters are the same as those of Fig. [3j The least-squares fit to these 
data yields v = 0.97(2) in the thermodynamic limit. A possible systematic error is considered in 
the text. 



Here, the operators {Qf} are given by the relations, Q] = (Sf)"^ — {SfY, Qf = [2(5*^)^ — 
{Sfy - {Sff]/VS, Qf = SfSf + SfSf, Qt = SfSt + SfSl and = S^Sf + S^Sf. This 
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reduction formula is a key ingredient in Sec. HIl 
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